Studies of the radiation environment on the Mars surface using the Geant4 toolkit* 
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The radiation environment on the surface of Mars is a potential threat for future manned exploration missions 
to this planet. In this study, a simple geometrical model was built for simulating the radiation environment on 
the Mars surface caused by galactic cosmic rays (GCRs); the model was built and studied using the Geant4 
toolkit. The simulation results were compared with the data reported by a radiation assessment detector (RAD). 
The simulated spectra of neutrons, photons, protons, œ particles, and particle groups Z = 3-5, Z = 6-8, 
Z = 9-13, and Z = 14-24 were in a reasonable agreement with the RAD data. However, for deuterons, 
tritons, and *He, the simulations yielded much smaller values than for the corresponding RAD data. In addition, 
the particles’ spectra within the 90° zenith angle were also obtained. Based on these spectra, we calculated the 
radiation dose that would have been received by an average human body on Mars. The distribution of the dose 
throughout the human body was not uniform. The absorbed and equivalent doses for the brain were the highest 
among all of the organs, reaching 62.0-+1.7 mGy/y and 234.1+8.0 mSv/y, respectively. The average absorbed 
and equivalent doses for the entire body were approximately 44 mGy/y and 153 mSv/y, respectively. Further 
analysis revealed that most of the radiation dose was owing to a particles, protons, and heavy ions. We then 
studied the shielding effect of the Mars soil with respect to the radiation. The body dose decreased significantly 
with increasing soil depth. At the depth of 1.5 m, the effective dose for the entire body was 17.9+2.4 mSv/y, 
lower than the dose limit for occupational exposure. At the depth of 3 m, the effective dose to the body was 
2.7+1.0 mSv/y, still higher than the accepted dose limit. 
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I. INTRODUCTION 


Space exploration has developed rapidly, and significant 
progress has been made in the past decades. Mars has be- 
come a hotspot for deep-space exploration missions. Mars 
exploration programs have been carried out, to learn about 
the evolution of Mars and for determining whether life ex- 
ists on Mars [1]. The manned mission to Mars has also been 
on schedule [2]. A key factor that must be considered for 
manned space exploration missions is the space radiation, 
which is caused by galactic cosmic rays (GCRs) and solar 
particle events (SPEs) [3-5]. In the space, energetic parti- 
cles can penetrate spacecraft and harm the instruments and 
astronauts [6]. Radiation-related health risks constitute a se- 
rious challenge for long-term manned space exploration mis- 
sions [7, 8]. GCRs are mainly composed of charged nu- 
clei, such as protons, a particles, and heavier nuclei with 
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a broad range of energies, from sub-MeV to hundreds of 
TeV [9, 10]. SPEs mainly consist of protons and a particles, 
with relatively lower energies compared with GCRs. Owing 
to the planet’s thin atmospheric layer and the lack of a global 
magnetic field [11], the radiation intensity on Mars is much 
higher than that on Earth. Knowing the radiation environ- 
ment on Mars and the radiation dose that is likely to be re- 
ceived by astronauts is very important for future manned mis- 
sions to Mars [11, 12]. This knowledge will allow decision- 
makers to accurately assess the corresponding health risks 
and is likely to provide a reference for designing radiation- 
protection gear [13, 14]. 


The radiation environment on Mars is complex and is de- 
fined by several factors. The shielding effect of the Mars at- 
mosphere is much weaker than that of the Earth atmosphere. 
As a result, energetic particles in GCRs and SEPs can pene- 
trate the Mars atmosphere, reaching the planet’s surface. Ow- 
ing to the irregularity and infrequency of SPEs, only GCRs 
were considered in this study. The interaction between GCRs 
and atmospheric molecules generates many secondary parti- 
cles, which further interact with the planet’s atmosphere, re- 
sulting in the complex radiation spectrum on the Mars sur- 
face. Albedo particles that are produced on the planet’s sur- 
face also contribute to the radiation environment. The three 
main components of the radiation environment on the Mars 
surface were considered in this study. The Curiosity rover 
has carried out a long-term scientific exploration on Mars 
since it successfully landed on the planet in 2012. At the 
same time, an onboard radiation assessment detector (RAD) 


has been collecting radiation-related information on the Mars 
surface, including the particles’ spectrum and the absorbed 
dose rate [15]. 

Some studies have shown that long-term exposure to the 
space radiation can harm humans [2]. Radiation can di- 
rectly and indirectly affect DNA molecules in the human 
body, causing cell death or mutation and increasing the risk 
of cancer [16—18] and other diseases. The National Aero- 
nautics and Space Administration (NASA) categorized the 
health risks owing to the space radiation into four groups: 1) 
degenerative tissue effects, 2) carcinogenesis, 3) acute and 
late central nervous system (CNS) effects, and 4) radiation- 
related syndromes [19]. Accurate calculation of the radia- 
tion dose received by astronauts is necessary for assessing 
the health risks associated with the space radiation. How- 
ever, it is difficult to directly calculate the particles’ spec- 
tra on the Mars surface and the associated radiation dose. 
Monte Carlo simulations can be used for addressing this 
problem. Some transport codes have been used for study- 
ing the radiation environment in the space and on Mars, 
such as HZETRN [20, 21], HETC-HEDS [22], FLUKA [23], 
PHITS [24], MCNP6 [25, 26] and Geant4 [27-30]. In 
Ref. [20-26], Slaba et al. studied the radiation environment 
on the Mars surface using these tools. Most of these studies 
have mainly focused on simulating the radiation environment 
on the Mars surface. Our work addresses two aspects of this 
research area. First, we used the GEQMD model to deal with 
the inelastic scattering of GCRs in the Mars atmosphere and 
soil, and for obtaining the particles’ spectra on the Mars sur- 
face. Then, based on the simulated particles’ spectra, we cal- 
culated the radiation dose received by a typical human body. 
To our knowledge, this is the first study attempting to estimate 
the radiation dosage to individual organs of a typical human 
on the Mars surface. The results and conclusions of this study 
can be used when planning future manned Mars exploration 
missions. 


II. THE RADIATION ENVIRONMENT ON MARS 
A. Model setup 


To reduce the simulation time, a simple geometrical model 
consisting of a cylinder (radius, 1 km) was built, to repre- 
sent the Martian environment. The cylinder was divided into 
12 layers to represent the Martian regolith and atmosphere. 
The first layer from the bottom corresponded to the Martian 
regolith (thickness, 3 m; density, 1.7 g/cm?) [21]. The sec- 
ond layer was 4-m-thick, representing the atmospheric condi- 
tions on the Mars surface. From the third layer to the twelfth 
layer, each layer was 1-m-thick, and the column densities of 
the different layers were 11.820, 5.915, 2.624, 1.069, 0.395, 
0.131, 0.035, 0.009, 0.003, and 0.001 g/cm?, respectively. 
The overall column density of these layers was approximately 
22 g/cm? [21, 31, 32], representing the Mars atmosphere. The 
detector was a cylinder (radius, 1 km; thickness, 3 mm) and 
was positioned! m above the soil, for detecting and register- 
ing the Mars radiation environment. The planet’s soil and 


atmosphere components were derived from Ref. [21] and are 
listed in Table 1. Primary GCRs were incident from the top 
of the geometrical structure, in the isotropic fashion. Once 
a particle passed through the detector, the information about 
the event was recorded, including the particle’s type, its ki- 
netic energy, and its angle of incidence. 


Table 1. Modeled soil and atmosphere components of Mars [21]. 
" Martian atmosphere 
Martian atmosphere CO2(95.7%), N (2.7%), Ar(1.6%) 

: : S102 (51.2%),Fe20; (9.3%), H20 (7.4% 
SAA a o ~ 
The radiation environment simulations were performed us- 

ing the Geant4-10.5.1 package (version10.5.1). This code is 
a track structure code that was developed at CERN for sim- 
ulating the passage of particles through matter [29, 30]. It 
provides many models that can be chosen by the user accord- 
ing to the specific application, and it has been widely used 
in high-energy physics, medical physics, and radiobiology. 
The physical parameters that were used for simulating the 
interaction between GCRs and the Mars medium are listed 
in Table 2. The Emstandard_opt3 physics list was used for 
describing the electromagnetic interactions of photons, lep- 
tons, hadrons, and ions, with energies in the 1 keV—10 PeV 
range [33]. The G4HadronPhysicsQGSP_BIC_HP physics 
list [34] was used for describing the inelastic scattering of 
nucleons and mesons with nuclei, for energies in the 0 eV- 
100 TeV range. In this physics list, the quark gluon string 
model (QGS) was used for describing high-energy hadronic 
interactions (12 GeV—100 TeV range). The binary cascade 
model (BIC) was used for dealing with nucleon-nucleus inter- 
actions at intermediate energies, and high-precision neutron 
models (HPs) were used for dealing with neutron-nucleus in- 
teractions involving energies below 20 MeV. For the inelastic 
interactions of nuclei, the G4QMD model that was validated 
in our previous work was used [35]. In this model, all nucle- 
ons were considered as participant particles, and each nucleon 
was treated as a Gaussian wave packet. Equations of the parti- 
cles’ motion were derived based on the time-dependent vari- 
ation principle [36]. A detailed description of these models 
can be found in the Geant4 Physics Reference Manual[36]. 


Table 2. The physics list 


Models Description 

Emstandard_opt3 Electromagnetic model 
G4HadronPhysicQGSP_BIC_HP Hadronic inelastic model 
G4IonQMDPhysics Hadronic inelastic model for ions 


In this study, only GCRs were considered as contributing 
to the radiation environment of Mars. GCRs mainly consist 
of protons, œ particles, and heavy ions. Although the pro- 
portion of heavy particles is very small, they are important 
for radiation protection, owing to their high relative biolog- 
ical effectiveness (RBE). In our simulations, particles with 
atomic numbers Z = 1-26 were considered for GCRs. The 
energy spectra of these particles were calculated using the 


GCR model developed by Matthia et al. [10]. In that model, 
a single parameter W was used for describing the modula- 
tion effect of solar activity [10]. A part of the primary GCR 
spectra used in the simulation is shown in Fig. 1. The mod- 
ulation parameter W was calculated according to the Oulu 
neutron monitor count rate (from http://cosmicrays.oulu.fi/) 
during the period from August 2012 to January 2013, and an 
average value (W = 56) was used in this work. 
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Fig. 1. (Color online) A part of the primary GCR spectra used in the 
simulation. 


B. Results 


When GCRs pass through a planet’s atmosphere, many 
neutrons and photons are induced owing to the interactions 
between the GCRs and atmospheric molecules. The simu- 
lated neutron and photon spectra on the Mars surface are plot- 
ted in Fig. 2, together with the experimental data [21, 37]. 
As shown in this figure, the simulation results for neutrons 
were slightly higher than the corresponding experimental 
data, for energies below 100 MeV/n, while the simulation 
results for photons agreed well with the RAD data. For neu- 
trons and photons, the average relative deviations of the simu- 
lation results from the corresponding experimental data were 
60.8% and 56.0%, respectively. The simulation results for 
the charged particles’ spectra, for zenith angles smaller than 
30°, are plotted in Figs. 3 and 4. Heavy particles can be di- 
vided into several groups, for proper comparison to the RAD 
data[15]. Evidently, the simulated spectra of protons, œ par- 
ticles, and particle groups Z = 6 — 9, Z = 9 — 13, and 
Z = 14 — 24 were consistent with the corresponding RAD 
data, and the corresponding average relative deviations were 
27.7%, 37.9%, 40.9%, 60.7%, and 40.1%. For the particle 
groups Z = 3— 5 and Z > 24, the simulation results were 
slightly lower than the measured data. However, the simu- 
lated spectra of deuterons, tritons, and 3He were much lower 
than the corresponding RAD data, and the average relative 
deviation was approximately 90%. These particles are pure 
secondary particles generated in reactions between primary 
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Fig. 2. (Color online) Simulated and measured photon and neutron 
spectra, for the Mars surface. 


GCRs and atmospheric molecules, and the observed discrep- 
ancy might be owing to the lower generation cross-sections. 
Based on the previous work, the cross-sections for light frag- 
ments simulated by the G4QMD model were lower than the 
experimental data[35]. 
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Fig. 3. (Color online) Simulated and measured charged particles’ 
spectra, for zenith angles smaller than 30°: (a) protons, deuterons, 
and tritons; (b) alpha particles and *He. 


Heavy charged particles incident from the entire upper 
hemisphere were also recorded in the simulation, and the sim- 
ulated spectra for zenith angles smaller than 90° are plotted 
in Fig. 5. Compared with the spectrum within the 30° zenith 
angle, for light particles, the flux averaged over the 90° zenith 
angle was close to the flux for zenith angles below 30°. How- 
ever, for heavy particles, the flux for the zenith angles below 
90° was slightly smaller than the flux for zenith angles be- 
low 30°. The spectra of electrons, muons, and pions were 
also calculated, owing to their significant contributions to the 
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Fig. 4. (Color online) Simulated and measured charged particles’ 


spectra, for zenith angles smaller than 30°: (a) Li-Be-B and C-N-O; 
(b) particle groups Z = 9-13, Z = 14-24, and Z > 24. 


absorbed dose, and are plotted in Fig. 6. 
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Fig. 5. (Color online) Simulated spectra of particles, for zenith an- 
gles smaller than 90°: (a) protons, deuterons, tritons, “He, and alpha 
particles; (b) Li-Be-B, C-N-O, F-Al, Si-Cr, and Mn-Fe. 


HI. RADIATION DOSE RECEIVED BY A TYPICAL 
HUMAN BODY 


Accurate estimation of the dose received by astronauts is an 
essential precondition for successful manned deep-space ex- 
ploration missions. In addition, with the wide use of nuclear 
technology in material physics and medical physics [38—46], 
it is also important to be able to calculate the radiation dose 
received by a typical human in the corresponding radiation 
environment. Based on the prior simulation results, we cal- 
culated the radiation dose that would have been received by 
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Fig. 6. (Color online) Simulated spectra of electrons, pions, and 
muonsm, for zenith angles smaller than 90°. 


astronauts on the Mars surface. 

According to the International Commission on Radiologi- 
cal Protection (ICRP) [47], the equivalent dose for a tissue, 
Hr, is defined by 


Hr =X _ WrDrr (1) 
R 


where Drp is the absorbed dose for tissue T caused by radia- 
tion R, and Wp is the radiation weighting factor for radiation 
R. The effective dose, Æ, is defined as 


E= So WrHr (2) 
T 


where Wr is the tissue weighting factor for tissue T. Table 
3 and Table 4 list the values of Wr and Wp, respectively. 
The ICRP recommends that the effective dose received dur- 
ing occupational exposure and in public should not exceed 20 
mSv/y and 1 mSv/y [47], respectively. 


Table 3. Tissue weighting factors Wr (extracted from Ref [47]). 
Wr Wr 


Tissue 
Bone-marrow (red), colon, 
Lung, stomach, breast, remainder tissues 0.12 
Gonads 0.08 
Bladder, oesophagus, liver, thyroid 0.04 
Bone surface, brain, salivary glands, skin 0.01 


0.72 
0.08 
0.16 
0.04 


A. Model setup 


The geometrical simulation model is shown in Fig. 7. The 
human phantom model [48] was borrowed from the Geant4 
advanced example “human phantom” and placed on the Mars 
surface. A cubical regolith (10x 106 m?) was built, with the 
density of 1.7 g/cm. Then, the simulated particles’ spectra 
within the 90° zenith angle (in Figs. 5—6) were used for sim- 
ulating primary particles and projecting them onto the surface 


Table 4. Radiation weighting factors Wp (extracted from ref [47]). 


Radiation type (R) 


Radiation weighting factors (WR) 


Photons 
Electrons and muons 
Protons and charged pions 
a particles, fission 
fragments, heavy ions 


1 
1 
2 


20 


2.5 + 18.2 exp {—[In(F,)]?/6}, En < 1 MeV 


Neutrons 


5.0 + 17.0 exp {—[In(2E;)]?/6}, 1 MeV < En < 50 MeV 


2.5 + 3.25 exp {—[In(0.04E,)]?/6}, Ea > 50 MeV 


from the top of the geometrical model. For each step in the 
simulation, the energy deposited by different particles into 
different human body organs was recorded. After the simu- 
lation, the values of the total deposited energy, the absorbed 
dose, and the equivalent dose delivered by different particles 
to different organs were calculated. In the simulation, the to- 
tal sampling number of each primary particle was set to 5 mil- 
lion. For each scenario, the simulation was repeated six times, 
and the uncertainties shown in this paper represent, in each 
case, one standard deviation from the corresponding mean. 


Fig. 7. (Color online) Geometrical simulation model. 


B. Results 


Figure 8 shows the radiation dose received by different hu- 
man body organs, for a typical human staying on the Mars 
surface for one year. The distribution of the deposited en- 
ergies in the body organs is plotted in Fig. 8(a). Ev- 
idently, the energy deposited in the trunk was the high- 
est, reaching (9.05 + 0.04) x 10t? MeV/y. The body part 
receiving the second highest amount of the radiation en- 
ergy were legs, with the deposited energy per leg reach- 
ing (2.02 + 0.02) x 10t? MeV/y. This large deposition oc- 
curred because the mass and volume of the trunk and legs 
were larger than those of the other body parts. The amount 
of energy deposited in the head and skull was also high; how- 
ever, the amounts of energy deposited in smaller organs, such 
as left/right adrenal and thyroid glands, wer among the low- 
est. Figure 8(b) shows the simulated absorbed dose distribu- 
tion across the different body parts. The distribution of the 
absorbed dose differs significantly from the distribution of 
deposited energy. The absorbed dose was the highest for the 


brain, reaching 62.0+1.7 mGy/y, followed by the thyroid and 
skull. The absorbed doses for the thymus, upper spine, left- 
/right lung, and left/right clavicle were also relatively high; 
however, the absorbed doses for the lower body organs were 
relatively low; for example, the absorbed doses for the left- 
/right leg and leg bone were approximately 35 mGy/y. Evi- 
dently, the difference between the doses received by the dif- 
ferent organs was significant, because primary particles were 
incident from above; consequently, some smaller-energy par- 
ticles only affected the superficial layers of the human body. 
Thus, upper-body organs, such as the brain and the thyroid 
gland, received higher radiation doses than lower-body or- 
gans. The average absorbed dose for the entire body was 
43.8+0.4 mGy/y. Taking the quality of different radiation 
sources into account, we calculated the equivalent doses for 
the different organs, and the results are shown in Fig. 8(c). 
The equivalent dose for the brain was the highest, reaching 
234.1+8.0 mSv/y, agreeing well with the corresponding RAD 
data [31]. For the skull, adrenal, and upper spine, the equiv- 
alent doses were also relatively high, all above 200 mSv/y. 
For the leg bone, the equivalent dose was the lowest, just 
116.4+1.7 mSv/y. The average equivalent dose for the en- 
tire body was approximately 153 mSv/y, much lower than the 
corresponding RAD data. We posit that this difference could 
be attributed to the fact that a typical human body is much 
larger than the RAD detector, and because the mass distribu- 
tion in the body is not uniform. The lower part of the body 
and inner organs received less radiation, so that the average 
equivalent dose for the entire body was lower than the corre- 
sponding RAD data. 


To analyze the contributions of different particles to the ra- 
diation dose received by the human body, we extracted the de- 
posited energy and the absorbed dose in the human body ow- 
ing to electrons, photons, pions, muons, protons, a particles, 
and heavy ions (Z > 3), respectively. The results are plotted 
in Fig. 9, except for photons, because their contribution was 
far less than those of other particles. The energy deposited by 
protons was (1.08-+0.01)x 1018 MeV/y, accounting for 56.1% 
of the total energy deposited in the human body. The second 
largest contribution was that of electrons; the corresponding 
deposited energy accounted for 22.6% of the total energy de- 
posited in the body. The energies depositions by a particles, 
heavy ions, pions, and muons were smaller, accounting for 
6.1%, 3.6%, 6.6%, and 3.9% of the total deposited energy, re- 
spectively. The contribution of neutrons to the total deposited 
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Fig. 8. (Color online) Radiation received by a human body on the Mars surface: (a) distribution of the radiation energy deposition in the body 
organs; (b) distribution of the absorbed dose in the body organs; (c) distribution of the equivalent dose in the body organs. The horizontal co- 
ordinate represents different organs: 1-Brain, 2-Head, 3-Heart, 4-LeftAdrenal, 5-LeftArmBone, 6-LeftClavicle, 7-LeftKidney, 8-LeftLeg, 9- 
LeftLegBone, 10-LeftLung, 11-LeftScapula, 12-LeftTeste, 13-Liver, 14-LowerLargeIntestine, 15-MaleGenitalia, 16-MiddleLowerSpine, 17- 
Pancreas, 18-Pelvis, 19-RibCage, 20-RightAdrenal, 21-RightArmBone, 22-RightClavicle, 23-RightKidney, 24-RightLeg, 25-RightLegBone, 
26-RightLung, 27-RightScapula, 28-RightTeste, 29-Skull, 30-SmallIntestine, 31-Spleen, 32-Stomach, 33-Thymus, 34-Thyroid, 35-Trunk, 


36-UpperLargelIntestine, 37-UpperSpine, 38-UrinaryBladder. 


energy was only 0.07%, much lower than that of protons. 
This was because the reaction cross-section of neutrons with 
matter was very small, and most of the deposited energy was 
owing to secondary particles. The contributions of the differ- 
ent particles to the absorbed dose were similar to those for the 
deposited energy. However, for the equivalent dose, the situ- 
ation was different. œ particles contributed mostly to the total 
equivalent dose, which reached 53.7+1.3 mSv/y, accounting 
for 35.1% of the total equivalent dose. The equivalent doses 
Owing to protons and heavy ions were 49.1+0.6 mSv/y and 
31.9+0.7 mSv/y, accounting for 32.1% and 20.9% of the total 
equivalent dose, respectively. The contributions of electrons, 
pions, muons, and neutrons were slightly smaller, accounting 
for only 6.5%, 3.8%, 1.4%, and 0.24% of the total equivalent 
dose, respectively. Although the deposited energies owing to 
q particles and heavy ions were not very high, they were the 
main contributors to the equivalent dose. This was because 
the radiation weighting factors of a particles and heavy ions 
were much larger than those of the other particles. Therefore, 


it is necessary to consider heavy ions in GCRs, when design- 
ing future manned deep-space exploration missions. 


It can be observed from Figs. 2-4 that the simulated spec- 
tra of deuterons, tritons, and 3He were much lower than those 
of the measured data. To assess the effect of the differences 
on the calculated doses to the human body, we multiplied the 
spectra by a scaling factor, for matching with the measured 
data, and recalculated the dose. The absorbed dose and the 
dose equivalents for the body organs, based on the scaled 
spectra, increased by approximately 22% on average, com- 
pared with the previous results. The average absorbed dose 
and the dose equivalent for the human body were 54.0 mGy/y 
and 191.6 mSv/y, respectively. The results show that these 
particles strongly affect the dose calculation, and the models 
in the Geant4 package should be improved for better describ- 
ing nucleus-nucleus interactions. 
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Fig. 9. (Color online) Contributions of different particles to the radiation dose received by the human body. 


IV. THE SHIELDING EFFECT OF THE MARTIAN SOIL 


The simulation results suggest that the dose that would 
have been received by astronauts on the Mars surface is much 
higher than that received on the Earth surface. Therefore, ap- 
propriate radiation protection measures are necessary. When 
astronauts will carry out long-term exploration missions on 
the Mars surface, it is likely that they will use Martian soil 
to build shelters, for protecting themselves from radiation. 
Thus, the shielding effect of Martian soil is worth investigat- 
ing. 


A. Model setup 


A hollow cube (5x5x3 m) was built and placed in the 
Mars soil, and the human phantom was placed in the cube. 
The thickness of the soil above the human phantom, denoted 
by d, was adjustable. In the simulation, we considered 9 lev- 
els of the soil thickness, namely d = 0.1 m, 0.3 m, 0.5 m, 
0.7 m, 1.0 m, 1.5 m, 2.0 m, 2.5 m, and 3 m; for each one of 
these scenarios we computed the dose received by the human 
body. The primary particles were projected from the top of 
the soil according to the simulated particles’ spectra on the 


Mars surface. The tissue energy deposition owing to different 
particles was recorded. 


B. Results 


The radiation doses on the human body at different depths 
under the Mars surface are plotted in Fig. 10. The deposited 
energy and dose rate on the Mars surface were also included 
for comparison. Evidently, with increasing the soil depth, the 
deposited energy, absorbed dose, and equivalent dose in the 
organs all decreased rapidly. At the depths of 0.1 m, 0.5 m, 
and 1.0 m, the absorbed dose was lower by 7.6%, 41.3%, and 
71.1%, respectively, compared with the dose on the Mars sur- 
face. The reduction rate of the absorbed dose decreased as 
the soil thickness continued to increase. As shown in Fig. 
10(c), the equivalent dose in the organs dropped faster than 
the absorbed dose at small depths. For example, at the depth 
of 0.1 m, the equivalent doses on the brain and whole body 
decreased by 30.6% and 28.1%, respectively. We attribute 
this to the fact that energetic heavy ions collided with the 
Mars soil and broke into lighter fragments, with smaller ra- 
diation weighting factors. The variation of the effective dose 
with the soil depth was similar to that of the dose equivalent 
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Fig. 10. (Color online) The dose on the human body, at different 
depths under the Mars surface. 


and at the soil depths of 0.1 m, 0.5 m, and 1.0 m, the ef- 
fective dose on the entire body was 121.0+9.7 mSv/y, 67.6 
+7.9 mSv/y, and 35.0+2.4 mSv/y, respectively. At the depth 
of 1.5 m, the average equivalent dose and effective dose on 
the entire body were 17.4+0.5 mSv/y and 17.9+2.4 mSv/y, 
respectively, lower than the dose limit of 20 mSv/y for occu- 
pational exposure. At the depth of 3 m, the equivalent doses 
received by different organs were between 2.0 mSv/y and 3.8 
mSv/y, and the effective dose on the entire body was 2.7+1.0 
mSv/y, still higher than the dose limit for public exposure. To 
confine the effective dose below 1 mSv/y, the thickness of the 
soil needs to be larger than 3 m. 


V. CONCLUSIONS 


In this study, we first simulated the radiation environment 
on the Mars surface using the Geant4 toolkit. In the simula- 
tion, a simple geometrical model and G4QMD physics model 
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